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Abstract. Matter bounces are bouncing scenarios wherein the universe contracts as in a 
matter dominated phase at early times. Such scenarios are known to lead to a scale invariant 
spectrum of tensor perturbations, just as de Sitter inflation does. In this work, we examine if 
the tensor bi-spectrum can discriminate between the inflationary and the bouncing scenarios. 
Using the Maldacena formalism, we analytically evaluate the tensor bi-spectrum in a matter 
bounce for an arbitrary triangular configuration of the wavevectors. We show that, over 
scales of cosmological interest, the non-Gaussianity parameter that characterizes the 
amplitude of the tensor bi-spectrum is quite small when compared to the corresponding 
values in de Sitter inflation. During inflation, the amplitude of the tensor perturbations 
freeze on super-Hubble scales, a behavior that results in the so-called consistency condition 
relating the tensor bi-spectrum and the power spectrum in the squeezed limit. In contrast, 
in the bouncing scenarios, the amplitude of the tensor perturbations grow strongly as one 
approaches the bounce, which suggests that the consistency condition will not be valid in 
such situations. We explicitly show that the consistency relation is indeed violated in the 
matter bounce. We discuss the implications of the results. 
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1 Introduction 

Bouncing models correspond to situations wherein the universe initially goes through a pe¬ 
riod of contraction until the scale factor reaches a certain minimum value before transiting 
to the expanding phase. They offer an alternative to inflation to overcome the horizon prob¬ 
lem, as they permit well motivated, Minkowski-like initial conditions to be imposed on the 
perturbations at early times during the contracting phase (see, for instance. Refs. [1-18]; 
for reviews, see Refs. [19-21]). Interestingly, certain bouncing scenarios can lead to nearly 
scale invariant perturbation spectra (see, for example. Refs. [3, 11, 18]), as is required by 
observations [22, 23]. For instance, a bouncing model wherein the universe goes through a 
contracting phase as driven by matter—referred to as a matter bounce—is known to lead 
to an exactly scale invariant spectrum of tensor perturbations as in de Sitter inflation [1- 
3]. Clearly, it will be worthwhile to examine if non-Gaussianities can help us discriminate 
between such scenarios [24-26]. 

The most dominant of the non-Gaussian signatures are the non-vanishing three-point 
functions involving the scalars as well as the tensors [27-34]. In order to drive a bounce, it is 
well known that one requires matter fields that violate the null energy condition. Therefore, 
analyzing the evolution of the scalar perturbations require suitable modelling of the matter 
fields [19-21]. In contrast, the tensor perturbations depend only on the scale factor and 
hence are simpler to study. For this reason, we shall focus on the tensor bi-spectrum in this 
work. Further, we shall assume a specific functional form for the scale factor and we shall 
not attempt to construct sources that can give rise to such a behavior. 

An interesting aspect of the three-point functions is their property in the so-called 
squeezed limit wherein the wavelength of one of the three modes involved is much larger 
than the other two [27, 35-43]. In this limit, under fairly general conditions, it is known that 
the three-point functions can be expressed completely in terms of the two-point functions, 
a relation that is referred to as the consistency condition. We should mention that, while 
the scalar consistency relation has drawn most of the attention, it has been established that 
all the four three-point functions involving scalars and tensors satisfy similar relations under 
certain conditions [28-30, 44, 45]. It is interesting to examine if the three-point functions 
generated in the bouncing scenarios satisfy the consistency condition. In the context of infla¬ 
tion, it is well known that the consistency relations arise due to the fact that the amplitude 
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of the long wavelength mode freezes on super-Hubble scales. In contrast, in a bouncing uni¬ 
verse it can be readily shown that the amplitude of the long wavelength mode grows sharply 
as one approaches the bounce during the contracting phase. This behavior suggests that 
the consistency relation may not hold in bouncing models [24]. The primordial consistency 
conditions lead to corresponding imprints on the anisotropies in the cosmic microwave back¬ 
ground (see, for instance. Refs. [46-48]; in particular, see Ref. [49] for the signatures of the 
tensor modes) and the large scale structure (see, for example. Refs. [50-52]). It is clear that 
the consistency condition, if it can be confirmed by the observations, can help us discriminate 
between models of the early universe. 

The most comprehensive formalism to study the generation of non-Gaussianities in the 
early universe is the approach due to Maldacena [27] . In this work, we analytically evaluate 
the tensor bi-spectrum in a matter bounce using the Maldacena formalism. To arrive at 
the tensor bi-spectrum analytically, one requires not only the behavior of the tensor modes, 
one also needs to be able to evaluate a certain integral involving the scale factor and the 
tensor modes. We conveniently divide the evolution into three domains and use the analytic 
solutions available in these domains to carry out the integrals and obtain the tensor bi¬ 
spectrum. 

This paper is organized as follows. In the following section, considering a specific form 
for the scale factor, we shall divide the period before the bounce into two domains, the first 
corresponding to early times during the contracting phase and the other close to the bounce. 
We shall describe the analytic solutions to the tensor modes during these two domains and 
also discuss the behavior of the modes after the bounce to arrive at the corresponding tensor 
power spectrum over wavenumbers much smaller than the wavenumber associated with the 
bounce. We shall also compare the analytical solutions for the tensor modes with the cor¬ 
responding results obtained numerically. In Sec. 3, we shall quickly summarize the essential 
expressions describing the tensor bi-spectrum in the Maldacena formalism. We shall also 
introduce the tensor non-Gaussianity parameter a dimensionless quantity that reflects 
the amplitude of the tensor bi-spectrum. In Sec. 4, we shall evaluate the tensor bi-spectrum 
using the analytic solutions to the modes and the behavior of the scale factor in the three 
domains. We shall calculate the bi-spectrum for an arbitrary triangular configuration of the 
wavevectors. In Sec. 5, we shall illustrate the results in the equilateral and the squeezed 
limits. We shall show that the non-Gaussianity parameter that characterizes the tensor 
bi-spectrum is very small for cosmological scales and that the consistency relation is violated 
in the squeezed limit. We shall conclude with a brief discussion in Sec. 6. In an Appendix, we 
shall briefly outline a proof of the consistency condition satisfied by the tensor bi-spectrum 
during inflation. 

Note that we shall work with natural units wherein h = c = 1, and define the Planck 
mass to be Mpj = (SttG)”^/^. 

2 The tensor modes and the power spectrum 

We shall consider the background to be the spatially flat Priedmann-Lemaitre-Robertson- 
Walker (FLRW) metric that is described by the line-element 

ds^ = a^(??) (—dr/^ -|- bij d®* d®-^) , (2.1) 

where a(r/) denotes the scale factor and r] is the conformal time coordinate. We shall assume 
that the scale factor describing the bouncing scenario is given in terms of the conformal time 
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coordinate rj by the relation 

a{v) = «o (1 + v'^/Vo) = oo (1 + /cq V^) , 

where oq is the minimum value of the scale factor at the bounce (i.e. when 77 = 0) and 
Vo = denotes the time scale that determines the duration of the bounce. Note that, 
at very early times, viz. when rj <C — 77 O) the scale factor behaves as in a matter dominated 
universe {i.e. as a oc 77 ^) and, for this reason, such a bouncing model is often referred to as 
the matter bounce. We shall assume that the scale associated with the bounce, viz. ko, is of 
the order of the Planck scale Mpj. Therefore, the wavenumbers of cosmological interest are 
50-60 orders of magnitude smaller than the wavenumber ko. 

Upon taking into account the tensor perturbations characterized by 7 ^-, the spatially 
flat FLRW metric can be expressed as [27] 


ds2 


a^iv) 





d®* dx^ 




(2.3) 


Recall that the primordial perturbations are generated due to quantum fluctuations. On 
quantization, the tensor perturbation 7 ^ can be decomposed in terms of the corresponding 
Fourier modes hk as follows: 

= E/ -0m + ■ ( 2 . 4 ) 


In this decomposition, the pair of operators (a|,, 0) represent the annihilation and creation 
operators corresponding to the tensor modes associated with the wavevector k, and they 
satisfy the standard commutation relations. The quantity e|j(fc) represents the polarization 
tensor of the gravitational waves with their helicity being denoted by the index s. The 
transverse and traceless nature of the gravitational waves leads to the conditions efj(fc) = 
ki£fj{k) = 0. We shall work with a normalization such that £lj{k) efj{k) = 2 [27]. The 

tensor power spectrum, viz. Pj,(k}, is defined as follows: 


{0iVe)lLiVe)) 


{2^f 

2fc3 4 


V.^{k) (5®(fc + fc'), 


(2.5) 


where the expectation values on the left hand sides are to be evaluated in the specified initial 
quantum state of the perturbations, and r]e denotes a suitably late conformal time when the 
power spectrum is to be evaluated. The quantity R^- is given by [28-30, 44, 45] 

= E 4(^)^-n(fc)- (2-6) 


The tensor spectral index n.^ is defined as 




dlnRp(A:) 

dln/c 


The tensor modes hk satisfy the differential equation 



(2.7) 


( 2 . 8 ) 
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where the overprimes denote differentiation with respect to the conformal time rj. If we write 
hk = Uk/a, then the modes Uk satisfy the equation 


^k + (k - — ] Uk = 0. 


(2.9) 


The quantity a"/a corresponding to the scale factor (2.2) is given by 


^ _ 2kl 

a l + k^rf' 


( 2 . 10 ) 


which has essentially a Lorentzian profile. Note that the quantity a”/a exhibits a maximum 
at the bounce, with the maximum value being of the order of k^. Also, it goes to zero as 
T] —>■ ±oo. For modes of cosmological interest, one finds that S> a"/a at suitably early 
times {i.e. as r/ —)• —oo). In this domain, the quantity Uk oscillates and we can impose the 
standard initial conditions on these modes and study their evolution thereafter. 

Let us divide the period before the bounce into two domains, a domain corresponding 
to early times and another closer to the bounce. Let the first domain be determined by the 
condition —oo < rj < —arjo, where a is a relatively large number, which we shall set to be, 
say, 10®. The second domain evidently corresponds to —arjo < rj < 0. In the first domain, 
we can assume that the scale factor behaves as 


a{r]) ~ aoklrj^, 


( 2 . 11 ) 


so that a"/a = 2/r/^. Since the condition k'^ = a"/a corresponds to, say, % = the 

initial conditions can be imposed when ij rjk- The modes hk can be easily obtained in such 
a case and the positive frequency modes that correspond to the vacuum state at early times 
are given by [1-3] 


hk 


Mpj aok'^T]'^ \ krjJ 


( 2 . 12 ) 


Let us now consider the behavior of the modes in the second domain, i.e. when —arjo < 
ij < 0. Since we are interested in scales much smaller than ko, we shall assume that rjk <C 
—arjQ, which corresponds to the condition k <C k^ja. Therefore, in this domain, for scales 
of cosmological interest, the equation governing the tensor mode hk reduces to 


hl + — h't^a. 

a 

This equation can be immediately integrated to yield 

a^iv* 


h'kiv) - h'kiv*) 


a?‘{rj) 


(2.13) 


(2.14) 


where r/* is a suitably chosen time and the scale factor a(r/) is given by the complete expres¬ 
sion (2.2). On further integration, we obtain that 


rv (Iff 

hkiv) =hk{v*)+h'f,{ri^)a‘^{rj^) J (2.15) 

where we have chosen the constant of integration to be /ifc(??*)• If we choose r/* = —arjQ, 
we can make use of the solution (2.12) to determine /ifc(T/*) and h'f,{rj^). Note that, in the 
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domain of interest, the first term in the above expression is, evidently, a constant, while the 
second term grows rapidly as one approaches the bounce. Upon using the form (2.2) of the 
scale factor, we find that we can express the behavior of the mode hk in the second domain 
as 


hk = Ak + Bkfikorj), 


(2.16) 


where 

l + k^T]^ 

while the quantities Ak and Bk are given by 


Ak 


Bk 


V2 1_ 1_ ^ 

Mpj ao V « ^ / 

1 1 h I 2^2 

Mp, V^2aoa2 ^ ’ 


e*“"/"“ + i?A:/(«), 

/3iko 3 ik\ 
\a‘^ k a ko J 


^iak/ko 


(2.17) 


(2.18) 

(2.19) 


Let us now turn to the third domain, i.e. immediately after the bounce. In this case 
too, for modes such that k <C feo/a, the solution to hk is given by Eq. (2.16). We should 
highlight the fact that, whereas the bounce (2.2) is a symmetric one, the solution (2.16) 
is asymmetric in r]. Moreover, one may have naively expected the amplitude of the long 
wavelength modes to freeze once the universe starts expanding. This is largely true, though 
not completely so, and the behavior can possibly be attributed to the specific form of the 
scale factor (2.2). Note that, during this domain, while the first term in f{kor]) decays, the 
second term actually grows, albeit extremely mildly. We shall assume that, after the bounce, 
the universe transits to the conventional radiation domination epoch at, say, rj = j3r]o, where 
we shall set (3 ~ 10^. We should mention that this choice is somewhat arbitrary and we 
shall discuss the dependence of the tensor power spectrum and the bi-spectrum on f3 in due 
course. 

In order to understand the extent of accuracy of the approximations involved, it would be 
worthwhile to compare the above analytical results for the mode hk with the corresponding 
numerical results. Clearly, given the scale factor (2.2), it is a matter of integrating the 
differential equation (2.8), along with the standard Bunch-Davies initial conditions, to arrive 
at the behavior of hk- The conformal time coordinate does not prove to be an efficient time 
variable for numerical integration, particularly when a large range in the scale factor needs to 
be covered. In the context of inflation, it is often the e-fold N, defined as a{N) = oq exp N, 
where = 0 is a suitable time at which the scale factor takes the value oq, that is utilized to 
integrate the equations of motion (see, for instance, Refs. [53-56]). Due to the exponential 
factor involved, a small range in e-folds covers a large range in time and scale factor. However, 
since e^ is a monotonically increasing function, while it is useful to describe expanding 
universes, e-folds are not helpful in characterizing bouncing scenarios. In order to characterize 
a bounce, it would be convenient to choose a variable that is negative during the contracting 
phase of the universe, zero at the bounce and positive during the expanding phase. We shall 
choose to perform the integration using a new variable A7, which we call the e-N-fold, in 
terms of which the scale factor is defined as a{M) = ooexp (A72/2) [57]. We shall assume 
that N is zero at the bounce, with negative values representing the phase prior to the bounce 
and positive values after. 
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In terms of the e-N-fold, the differential equation (2.8) governing the evolution of the 
tensor modes can be expressed as 


CW2 + 


3M + 


1 dH 

H m 



d/tfc 

cW 



hk 0 , 


( 2 . 20 ) 


where H = a'/a'^ is the Hubble parameter. Given the scale factor (2.2), the corresponding 
Hubble parameter can be easily evaluated in terms of the conformal time rj. In order to 
express the Hubble parameter H in terms of the e-N-fold, we shall require rj as a function 
of M. Upon using the definition of the e-N-folds and the expression (2.2) for the scale factor, 
we obtain that 

r]{M) = ± - l) . (2.21) 

Since the Hubble parameter is negative during the contracting phase and positive during 
the expanding regime, we have to choose the root of ri{M) accordingly during each phase. 
From the expression for the Hubble parameter, we evaluate the coefficients of the differential 
equation (2.20) in terms of N. With the coefficients at hand, we numerically integrate 
the differential equation using a fifth order Runge-Kutta algorithm. We should mention 
that we have also independently checked the numerical results using Mathematica. Recall 
that, the initial conditions are imposed in a domain during the contracting phase wherein 
> a"/a. As is done in the context of inflation, we shall impose the initial conditions 
when = 10^ {o-''/o) corresponding to, say, the e-N-fold A/). It should be pointed out that 
the initial conditions on the different modes are imposed at different times. In terms of the 
e-N-folds, the standard Bunch-Davies initial conditions can be expressed as 


a(A/’j)\/^’ 

dhk iMi 

^ ~ a‘^{Mi) H{Mi) 

We impose these initial conditions well before the bounce and evolve the modes until a 
suitable time after the bounce. The tensor mode evaluated in such a manner has been 
plotted for a given wavenumber (such that fe/fco 1) in Fig. 1. The figure also contains a 
plot of the analytical results (2.12) and (2.16) for the same wavenumber. As is evident from 
the figure, prior to the bounce and immediately after, the analytical results match the exact 
numerical results exceedingly well. 

The tensor power spectrum after the bounce can be calculated using the solutions we 
have obtained. Recall that the tensor power spectrum is defined as 

'PT{k) = ^^\hk{r,)\\ (2.23) 

with the spectrum to be evaluated at a suitable time. If we evaluate the tensor power 
spectrum at r] = (3 tjq, we find that it can be expressed as 

rAk) = A^\A, + B,f{/3)\^ (2.24) 

Note that, a is a quantity that we have artificially introduced and the actual problem does 
not contain a. For A: <C ko/a and a sufficiently large a (as we had said, for a = 10® or so). 




a{Mi) 


(2.22a) 

(2.22b) 
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Figure 1. A comparison of the numerical results (in blue) with the analytical results (in red) for 
the amplitude of the tensor mode \hk\ corresponding to the wavenumber k/ko = 10 “^°. We have 
set ao = 10® and, for plotting the analytical results, we have chosen a = 10®. We have plotted the 
results from the initial e-N-fold A/) [when k^ = 10'^ (a"/a)] corresponding to the mode. While we have 
illustrated the exact numerical result till rather late times, we have plotted the analytical results until 
a time after the bounce when the power spectrum is evaluated (see discussion below). Evidently, the 
analytical and numerical results match extremely well, suggesting that the analytical approximation 
for the modes works to a very good accuracy. 


the above power spectrum reduces to a scale invariant form with a weak dependence on /?, if 
P is reasonably larger than unity. If we further assume that /3 is large enough (say, 10^), then 
the scale invariant amplitude is found to be: Vrj,{k) ~ 9/cg/(2Mpj Oq), as expected [1-3]. In 
Fig. 2, we have plotted the complete tensor power spectrum described by the expression (2.24) 
for a given set of parameters. We should stress that the power spectrum is actually valid 
only for modes which satisfy the condition k ko/a. It is evident from the figure that 
the power spectrum is strictly scale invariant over this domain. Moreover, we find that the 
spectrum indeed reduces to the above-mentioned scale invariant amplitude for small values 
of the wavenumbers. We have also evaluated the tensor power spectrum numerically using 
the method described above. We have computed the spectrum at a given time soon after the 
bounce (corresponding to /3 = 10^) for all the modes. We find that, for wavenumbers such 
that A: <C feo, the numerical analysis also leads to a scale invariant spectrum whose amplitude 
matches the above analytical result to about 1%. 
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Figure 2. The behavior of the tensor power spectrum has been plotted as a function of fc/fco for 
a wide range of wavenumbers. In plotting this figure, we have set fco/Mp, = 1, oq = 10^, a = 10® 
and /3 = 10^. We should emphasize that the approximations we have worked with are valid only over 
the domain wherein k <C ko/a. Clearly, the power spectrum is scale invariant in this domain. We 
also find that, at small wavenumbers, the tensor power spectrum has the expected scale invariant 
amplitude of Pj-ik) = 4.5 x 10“^° corresponding to ko/M^^ = 1 and oq = 10®. 

3 The tensor bi-spectrum and the corresponding non-Gaussianity 
parameter 

As we have mentioned, the most comprehensive formalism to calculate the three-point func¬ 
tions generated in the early universe is the formalism due to Maldacena [27]. The primary 
aim of Maldacena’s approach is to obtain the cubic order action that governs the scalar and 
the tensor perturbations using the ADM formalism. Then, based on the action, one arrives 
at the corresponding three-point functions using the standard rules of perturbative quantum 
field theory. 

The tensor bi-spectrum in Fourier space, viz. {ki, k 2 , k^), evaluated at 

the conformal time, say, r/e, is defined as 

^ (fci + h). 

(3.1) 

Note that the delta function on the right hand side implies that the wavevectors ki, k 2 and 
k^ form the edges of a triangle. For convenience, hereafter, we shall set 

5™™™(fci,fc2,fc3) = (2 7r)-9/2 (3.2) 

The tensor bi-spectrum Gl^^i^^'^'^^'^^^'^^{ki,k 2 ,k 3 ), calculated in the perturbative vacuum 
using the Maldacena formalism, can be written in terms of the modes as follows [27, 31- 
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34, 44]: 




/n^i ,.n''2 . n 

V mini,ij m 2 n 2 ,im 


jk2 


kz 

mznz,lj 


- 2 ^^n 2 ,ml ^mznz,ij) ^im hi + five permutations 

X [hk^{r]e) hk^ive) hkz{r]e)Gjjy{ki,k2,k3) 

+ complex conjugate], (3-3) 


where G-y-y-y{ki, k 2 , k^) is described by the integral 


g^^^{ki,k2,kh = -- / drj a'^hl^hl^hl^, 
^ Jr)i 


(3.4) 


with r][ denoting the time when the initial conditions are imposed on the perturbations and r/e 
representing the time when the bi-spectrum is to be evaluated. Also, we should mention that 
{hi-,k 2 i-,hi) denote the components of the three wavevectors {ki,k 2 ,k 3 ) along the i-spatial 
direction^. 

The dimensionless non-Gaussianity parameter that characterizes the amplitude of the 
tensor bi-spectrum is defined as [34] 

2 


Ki.{ki,k2,kh = - 


2h 


[kf kl ki fcs)] 


Trfcl TTfc2 J,3 

■‘■■‘■mini,mans ■‘■‘■m2n2,mn 


7^2 


1 -1 


tnn ^3 ’Pt(^ 2 ) + five permutations , (3.5) 


where the overbars on the indices imply that they need to be summed over all allowed values. 
Our aim in this work is to evaluate the magnitude and shape of the tensor bi-spectrum and 
the corresponding non-Gaussianity parameter and compare with, say, the results in de Sitter 
inflation. Therefore, for simplicity, we shall set the polarization tensor to unity. In such a 
case, the expression (3.3) for the tensor bi-spectrum above simplifies to 


^777(^11^2)^3) — “^pi [^fci (%) ^^2 (l?e) (l?e) ^777 (^1) ^2) ^3) 

-|- complex conjugate], (3-6) 

where the quantity Q^^^{ki, ^ 2 , h) is described by the integral 

i 

^ 777 (^ 1 , ^ 2 , ^ 3 ) = —T (^1 + ^2 + ^ 3 ) / dr/• (3-7) 

We shall choose r/j to be an early time during the contracting phase when the initial conditions 
are imposed on the modes (ie. when h ;§> a"/a) ^ r/e to be a suitably late time, say, some 

time after the bounce, when the bi-spectrum is evaluated. If we ignore the factors involving 
the polarization tensor, the non-Gaussianity parameter reduces to 


hj,hh,k2,k3) 



{ki,k2,k3)] 


2kl Vrj,{ki)Vrj,{k 2 ) + two permutations 


-1 


(3.8) 


^Such a clarification seems necessary to avoid confusion bet-ween fei, kz and kz which denote the wavenum¬ 
bers associated with the wavevectors fei, k 2 and kz, and the quantity ki which represents the component of 
the wavevector k along the r-spatial direction. 
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4 Evaluating the tensor bi-spectrum 


With the forms of the scale factor and the mode functions at hand, in order to arrive at the 
tensor bi-spectrum, it is now a matter of evaluating the integral (3.7) in the three domains. 

Let us begin by considering the first domain. Upon using the behavior (2.11) of the scale 
factor and the mode (2.12) in the first domain, we hnd that the quantity k 2 , k^) can 

be expressed as 


^ 777 (^ 1 ; ^ 2 , fcs) 


—i (/cf + k2 + A;|) 


4 ao kl y/ki /c2 fcs _ 


hikrj,,ko, a) + i 


1 1 1 

ki^ k2^ ks 


1 ^0 1 ^) 


1 




1 




1 


kl k2 k2 ks kl ks 


h{k^,ko,a) - 


kl k2 ^3 




, (4.1) 


where k^ = ki + k 2 + k^ and the quantities In{k^,ko,a) are described by the integrals 


In {k^ ^0; ^) 



r] 


(4.2) 


For n > 1, these integrals can be evaluated to yield 


Iri-\-l (^x 5 ^0? ^) 


1 

n 



^-iak^/ko 




i kr. 


n 


Ini^k^Y , ko , u) 


(4.3) 


while Ii{krj,, ko, a) is given by 


Ii{krj,, fco, a) = in + Ei(—ia krj^/ko), 


(4.4) 


where Ei(x) is the exponential integral function [58]. 

Let us now turn to evaluating G-y-y-y{ki,k 2 ,k 3 ) in the second domain. Upon using the 
behavior (2.2) of the scale factor and the mode (2.16), we find that the quantity can be 
expressed as 


g^^^{ki,k2,k 


3) = 


zoq [kl + k2 + fe|) 


4kn 




+ {Al Al Bl^ + Al^ Bl^ Al^ + Bl A*,^ A*J Ji(a) 

+ (Al Bl 54 + Bl Al 54 + Bl Bl^ Al^) J2{a) 

+BI, Bl^Bl^Ua) 


(4.5) 


where Jn{oi) are described by the integrals 

Jn{a) = [ da; [l + P{x), (4.6) 

J—a 

with the function f{x) being given by Eq. (2.17). The integrals Jo(a) and Ji(a) can be 
readily evaluated to obtain that 

T / \ ‘lo? ^ 

Jo(a) = a + ^- + — (4.7) 

3 5 


- 10 - 














and 




2\2 


n 2 (l + a^) ^ (l + a2) 


60 


15 


20 


8 a _i 
— tan a 
15 


— ^ “ T ^ a + :^ In (l + a^) . (4.8) 

15 5 15 

In contrast, the integrals J 2 {oi) and J'i{a) are more involved. The integral J 2 {ot) can be 
divided into three parts and written as 


J2(a) — J2l(a) + J22(a) + </23(«); 

where the integrals J 2 i(a) and J 22 (a) can be easily evaluated to be 

3 ’ 


>^ 21 ( 0 ) = / dxx^ = 

J —OL 

fO 

J 22 (a) = 2 / dxx (l + x^) tan 

J —OL 


(4.9) 


(4.10) 


2 I -t a \ _ 1 1 / _ 1 \ a 

= a 1 H-tan a - la — tan a) -. 

' 2 / 2 '6 


The quantity J23(a) is given by 

fO 2 , ' 

>^ 23 ( 0 :)= / dx (l + x^) (tan“^x)" 

J —a 

and, upon setting tan“^ x = y, it reduces to 


<^23(«) = 


dy y sec° y. 


' — tan ^ a 

The integral involved can be evaluated to be (see, for instance, Ref. [58]) 

/■ 2 6 -2/ (cosy - 2 y siny) 4 y (cosy - y siny) /II 8y^ 

/ dy y sec° y = -77^-^-7^-5-h ( — + 


(4.11) 


(4.12) 


(4.13) 


30 


10 cos® y 
15 


15 cos® y 


30 15 


tan y 


tan® y 16 ^ (- 1 )*^ 2^^ ( 2 ^" - l) 

-- ^ ~ --D 2 n, 


n=l 


(2n + 1) (2n)! 


(4.14) 


where B 2 n are the Bernoulli numbers. Needless to add, this result can be used to arrive at 
J 23 (a). We should add that the infinite series in the above expression is convergent, and we 
find that it can be expressed as follows [59]: 
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where C'(s, a) denotes the derivative of the Hurwitz zeta function ((s, a) with respect to the 
first argument and r(n) is the Gamma function. 

Let us now evaluate the last of the integrals, viz. Jsia). It proves to be convenient to 
divide the integral into four parts as follows; 


«^3(«) — <^3i(«) + '/ 32 (a) + J33(a) + 


(4.16) 


If we set tan ^ x = y, we find that the integrals J3 i(q;), Jziia) and J 33 {a) can be easily 
evaluated to be 


fO q,2 I 

'/31(a) = / dytan^y = ——+ -ln(l + a 2 ), 

/-tan-la / ^ 

/■O 


(4.17) 
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(4.18) 


'/33(a) — 3 


fO 


J — tan ^ a 

The integral J34(a) is given by 


dyy tan^ y sec^ y = — -a^ tan ^ a — - In (l + . (4.19) 
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J — tan“^ a. 


(4.20) 


' — tan ^ Q 

which can be evaluated using the result [58] 
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The infinite series in the above expression is convergent, and it can be expressed as [59] 
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(4.22) 


where, as before, C(s) is the Riemann zeta function, C,'{s,a) denotes the derivative of the 
Hurwitz zeta function C,{s,a) with respect to the first argument and r(n) is the Gamma 
function. 

Let us now consider the quantity Q^^^{ki,k 2 ,k^) in the third domain, i.e. from the 
bounce at r/ = 0 to ?] = /Sryo- In this domain, the modes hk and the scale factor have 
the same form as in the second domain. Therefore, it should be clear that, the quantity 
0^^^{ki, k 2 , fcs) too will be given by the expression (4.5), but with the integrals Jn{oi) being 
replaced by —Jni—P). 


5 Results 

We can now make use of the behavior of the mode hk at rje = j3rjo and substitute the 
results we have obtained above in the expressions (3.6) and (3.8) to arrive at the tensor bi¬ 
spectrum and the corresponding non-Gaussianity parameter for an arbitrary triangular 
configuration of the wavenumbers involved. The resulting expressions prove to be rather long 
and, for this reason, we shall illustrate the various results graphically for a set of suitable 
values of the parameters. Let us first compare the contributions from the three domains. 
Restricting ourselves to the equilateral limit, in Fig. 3, we have plotted the contributions to 
from the three domains that we have considered. It is evident from the figure that the 
contribution due to the third domain to the parameter turns out to be the maximum. 
We find that the third domain contributes the maximum in the squeezed limit as well. 

It is now a matter of adding the contributions from the three domains to arrive at 
the complete tensor bi-spectrum and the non-Gaussianity parameter In Fig. 4, we have 
plotted the behavior of the non-Gaussianity parameter in the equilateral and the squeezed 
limits. Three points concerning the figure require emphasis. To begin with, we should 
mention that the non-Gaussianity parameter behaves as k'^ in both the equilateral and 
the squeezed limits, with virtually the same amplitude. Secondly, the value of the parameter 
Is very small when compared to the values that occur in, say, de Sitter inflation wherein 
3/8 < ~ 1/2 (in this context, see Ref. [34]). Thirdly, since the tensor power spectrum 
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Figure 3. The contributions to the non-Gaussianity parameter in the equilateral limit from the 
first (in green), the second (in red) and the third (in blue) domains have been plotted as a function 
of fc/fco for a wide range of wavenumbers such that k kg/a. We have worked with the same set of 
values as in the previous figure. Clearly, the third domain gives rise to the maximum contribution to 
the non-Gaussianity parameter 


is strictly scale invariant for wavenumbers such that k <C koja, the amplitude of the non- 
Gaussianity parameter in the squeezed limit over such domain should be equal to 3/8, if 
the consistency relation holds true (see Appendix, also see Ref. [45]). Whereas, we find that 
is considerably smaller than 3/8 in the squeezed limit, which unambiguously implies that 
the consistency condition is violated. Evidently, this behavior can be attributed to the fact 
that the amplitude of the tensor mode does not freeze to a constant value at late times. 

At this stage, we need to discuss the dependence of the tensor power and bi-spectra on 
the parameters a and /3 that we have introduced. We find that the tensor power spectrum 
and the bi-spectrum do not significantly depend on a over a wide range of values, say, 
10'" < (X < 10^^. We had mentioned earlier that the tensor power spectrum has a rather 
weak dependence on /3. Whereas, we find that the tensor bi-spectrum grows roughly as 
However, /3 cannot be allowed to be too large for two reasons. One may a priori expect that 
the analytical approximation (2.16) will remain valid until the time —rjk = V^jk after the 
bounce. We had pointed out that the evolution of the mode is asymmetric in tj. Actually, 
it can be shown that (using numerical analysis) the analytical approximation (2.16) breaks 
down much before —rjk- For this reason, we have chosen (5 to be smaller than a. Moreover, 
a transition to the radiation dominated phase is expected to take place sometime after the 
bounce. It seems reasonable to expect that such a transition will occur when the scale factor 
is a ~ lO'^ oq, which corresponds to /3 = 10^. We find that our main conclusions, viz. that the 
value of is small over cosmological scales and that the consistency relation is violated in 
the squeezed limit, continue to remain valid even if we increase (5 by, say, a couple of orders 
of magnitude. 
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Figure 4. The behavior of the non-Gaussianity parameter in tti® equilateral (in blue) and the 
squeezed (in red) limits have been plotted as a function of fc/fco for a wide range of wavenumbers such 
that k <C ko/a. We have worked with the same set of values as in the earlier two figures. Clearly, 
the resulting is considerably small when compared to the values that arise in de Sitter inflation 
wherein 3/8 < < 1/2. Moreover, we find that behaves as k^ in the equilateral and the 

squeezed limits, with similar amplitudes. The fact that is much smaller than 3/8 in the squeezed 
limit implies that the consistency condition is violated. 


6 Discussion 

In this paper, we have analytically calculated the tensor bi-spectrum in a matter bounce 
using the Maldacena formalism. While the matter bounce leads to a scale invariant tensor 
power spectrum for scales of cosmological interest as de Sitter inflation does, we have shown 
that the non-Gaussianity parameter that characterizes the amplitude of the tensor bi¬ 
spectrum is much smaller than the corresponding values in de Sitter inflation. We have 
also shown that, due to the growth in amplitude of the tensor modes as one approaches the 
bounce, the consistency condition is not satisfied by the tensor bi-spectrum in the squeezed 
limit. Recall that, in the absence of detailed modelling of the bounce, we had assumed that 
ko ~ Mpi- We should however clarify that, since k ko for cosmological scales, our essential 
conclusions, viz. that is small and that the consistency condition is violated over such 
scales, will remain unaffected even if we choose ko to be a few orders of magnitude below 
the Planck scale. It will clearly be worthwhile to investigate these issues using a numerical 
approach in a wider class of bouncing models. 

In the bouncing scenario that we have considered, at very early times, i.e. during the 
first domain of our interest, the contribution to the non-Gaussianity parameter can be 
said to be small because the amplitude of the tensor perturbations themselves are small. In 
the second domain, although the scale factor decreases gradually to reach its minimum at 
the bounce, the non-Gaussianities become larger as the perturbations grow. In the third 
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domain, i.e. after the bounce, the scale factor increases steadily. Also, the amplitude of the 
perturbations do not freeze but grow slowly. Due to these reasons the contribution to the 
non-Gaussianity parameter is the largest from this regime. However, essentially due to the 
form of the scale factor, one finds that the parameter has an overall dependence. 
Since the scales of cosmological interest are about 50 to 60 orders below the Planck scale, 
the non-Gaussianity parameter proves to be very small over such scales. 

We believe that the results we have obtained have tremendous implications for the other 
three-point functions and, importantly, the scalar bi-spectrum. It seems clear that, due to the 
growth during the contracting phase near the bounce, the consistency conditions governing 
the other three-point functions will be violated as well [24]. This possibly can act as a 
powerful discriminator between the inflationary and bouncing scenarios. Within inflation, 
one requires peculiar situations to violate the consistency conditions [60, 61]. In contrast, in 
a bouncing scenario, the consistency relations seem to be violated rather naturally. Notably, 
situations involving violations of the consistency conditions have been considered as possible 
sources of spherical asymmetry in the early universe [28-30]. These aspects seem worth 
exploring in greater detail. 
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Appendix: The squeezed limit and the consistency condition 

An important property of the three-point functions is their behavior in the so-called squeezed 
limit [27, 35-43]. As we have discussed before, the squeezed limit corresponds to the situation 
wherein one of the three wavenumbers is much smaller than the other two. In such a limit, 
under certain conditions, it is known that all the three-point functions involving the scalars 
and tensors generated during inflation can be expressed entirely in terms of the two-point 
functions [28-30, 44, 45]. In the context of inflation, these consistency conditions arise 
essentially because of the fact that the amplitude of the long wavelength scalar and tensor 
modes freeze on super-Hubble scales. In this appendix, we shall outline the proof of the 
consistency condition satisfied by the tensor bi-spectrum during inflation. 

Since the amplitude of a long wavelength mode freezes on super-Hubble scales during 
inflation, such modes can be treated as a background as far as the smaller wavelength modes 
are concerned. Let us denote the constant amplitude of the long wavelength tensor mode 
as 7?-. In the presence of such a long wavelength mode, the background FLRW metric can 
be written as 

ds^ = —dt^ -|- a^(t) [e'*''^ ]ij d®* dx^, (.1) 

i.e. the spatial coordinates are modified according to a spatial transformation of the form 
x' = Ax, where the components of the matrix A are given by Ajj = [e'*' Under such a 

spatial transformation, the small wavelength tensor perturbation transforms as [45] 

7‘ ^det(A->)7y‘t (.2) 
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where det (A = 1. Under these conditions, we also obtain that 


IA ^ fc| = [1 - 7®. fcj fcj/(2fc^)] fc, (.3) 

where ki is the component of the wavevector k along the i-spatial direction and we have 
restricted ourselves to the leading order in 7?.. Moreover, one can show that 

5(3)(a-i + A-ifc2) = dQi{K)5^^\ki + k2) = 5^^\ki + k2), (.4) 


since det (A) = 1. Upon using the above results, we find that the tensor two-point function 
in the presence of a long wavelength mode can be written as 


(a 


fcl 7^2 \ 

m\ni Jm2n2/^ 


2kf 


mini ,m2n2 


V^{ki)6^^\ki + k2) 


1 - 


Urj, — 3 


% nii nij 


(.5) 


where hu = ku/ki and the long wavelength mode is denoted by the wavenumber k, while Urj, 
represents the tensor spectral index. The corresponding expression for the tensor bi-spectrum 
can be obtained from the above result to be 
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(. 6 ) 


where k^ has been considered to be the squeezed mode. The above relation wherein the 
tensor bi-spectrum has been expressed completely in terms of the power spectrum is known 
as the consistency condition [44, 45]. Upon substituting this expression in the definition for 
the tensor non-Gaussianity parameter [cf- Eq- (3-5)], we find that we can express the 
consistency relation in the squeezed limit as follows: 
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(•7) 


Actually, an overall minus sign occurs in the expression for in the squeezed limit due to 
the polarization tensors [45], Therefore, if we ignore the polarization tensors, in the domain 
where the tensor power spectrum is strictly scale invariant (i.e. when =0), the value of 
in the squeezed limit reduces to 3/8, if the consistency relation is satisfied. 
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